Stability of the Ni sites across the pressure-induced metanization in YNiO^ 
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The local environment of nickel atoms in YNiOs across the pressure- induced insulator to metal 
(IM) transition was studied using X-ray absorption spectroscopy (XAS) supported by ab initio 
calculations. The monotonic contraction of the NiOe units under applied pressure observed up 
to 13 GPa, stops in a limited pressure domain around 14 GPa, before resuming above 16 GPa. 
In this narrow pressure range, crystallographic modifications basically occur in the medium/long 
range, not in the NiOe octahedron, whereas the evolution of the near-edge XAS features can be 
associated to metallization. Ab initio calculations show that these features are related to medium 
range order, provided that the Ni-O-Ni angle enables a proper overlap of the Ni Cg and O 2p orbitals. 
Metallization is then not directly related to modifications in the average local geometry of the NiOa 
units but more likely to an inter-octahedra rearrangement. These outcomes provides evidences of 
the bandwidth driven nature of the IM transition. 

PACS numbers: 71.30.+h, 78. 70. Dm, 71.15.-m, 71.27.+a 



I. INTRODUCTION 

The Ni perovskite family RNiO^ {R = Y and rare-earth 
=/= La) presents a localized 3d-electron behavior at low 
temperatures and a first-order insulator to metal (IM) 
transition as temperature increaseJi^. The transition 
temperature (Tim) is determined by the R ionic radius 
that modifies the mismatch between the Ni-0 and R-0 
bond lengths. Such modification goes along with changes 
in the Ni-O-Ni superexchange angle. A decrease in Tjm 
reflects an increase in the Ni 3d eg and O 2p orbitals hy- 
bridization, concomitant with an increase in the band- 
widtlP'. In compounds with small R ions (Y, Ho, Er, Tm, 
Yb, and Lu) the crystallographic symmetry changes from 
monoclinic to orthorhombic across the IM transitiorP^^. 
The single Ni^~^ site in the metallic phase split in the in- 
sulating phase into two nonequivalent Nil and Ni2 sites, 
with slightly different average Ni-0 distances. These two 
different average distances are interpreted as a signature 
of charge order (or charge disproportionationpl^. It 
has also been proposed that the difference in Ni local 
environment reflects the presence of an ionic bonding at 
Jahn-Teller (JT) distorted Ni2 (larger) sites and cova- 
lent bonding at Nil (smaller) sites, with possible dy- 
namical fluctuations between both siteai^-li. Conflrm- 
ing such fluctuations, a strong electron-lattice coupling 
in a JT-distorted lattice has been evidenced in small- 
R compoundJ^. For the largest-i? compounds the time 
scale of the fluctuations is expected to be shorter and 
the splitting is hardly observable by elastic scattering 



technique j^ * ^° l, leading to the description of an average 
structure. Changes in the local symmetry can be tracked 
by X-ray absorption spectroscopy (XAS). 

As a structural technique able to probe dynamical 
changes, XAS has shown that the two Ni sites coexist 
in both in sulati ng and metallic state in several RNiOz 
compounds--3Ji^. In PrNiO^ local symmetry changes ac- 
companying electronic delocalization were found incom- 
patible with the average long range order proposed from 
X-ray and neutron diffractiorfl^. These studies empha- 
size the existence of an inhomogeneous structure at the 
local scale and suggest a common behavior in all RNiOs 
compounds for the local electronic and magnetic statdJ^l. 

Among the RNiO^ compounds, YNiO^ presents one of 
the largest monoclinic distortion at room temperature, 
due to the small size of the Y^^ ion^. The thermal- 
induced metallic phase occurs simultaneously with the 
vanishing of the long-range monoclinic distortiorP. How- 
ever, the two Ni sites coexistence, clearly observed by 
XAS even in the orthorhombic phase^'^ supports a 
model of stable short-range scale distortion. 

The application of an external pressure provides a unique 
tool to further investigate the relationship between struc- 
tural distortions and electronic properties. Hydrostatic 
pressure reduces the unit cell volume and shrinks Ni- 
O bond lengths while straig htening the bond angle and 
stabihzing the metaUic phasePE^lMl, in YNiO:^ Garcia 
Munoz and coworkers reported by X-ray diffraction a 
sudden structural modification around 14 GPa, consis- 
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tent with a monoclinic to orthorhombic transitiorP^. At 
this pressure, an increase in electronic conductivity and 
a phonon screening are also observed by infrared spec- 
troscopy. However, the limited resolution of the diffrac- 
tion experiments did not allowed a thorough study of 
the monoclinic distortion and the metallic phase was not 
confirmed. 

In the present paper we report in situ high pressure XAS 
experiments for YNiO^ up to 19 GPa. The expected IM 
transition is finger-printed by the XAS near-edge features 
around 14 GPa. Around that pressure we do not observe 
any modification in the NiOe geometry. The changes 
are mostly coming from a rearrangement of the octahe- 
dra that leads to a straightening of the superexchange 
angle. The subsequent conclusion is that the delocaliza- 
tion of the Cg electrons is due to band effects, enabled in 
the orthorhombic phase by favorable orbital alignments. 
Our results show that XAS spectroscopy may provide the 
essential information about octahedral links, hindered 
in scattering techniques by dynamical fluctuation in the 
short range scale. We show that the occurrence of the 
MI transition is not related to the exact local geome- 
try of the NiOg octahedra or charge disproportionation. 
It essentially depends on the middle range organization. 
This confirms the paramount importance of the octahe- 
dral tilting in the physics of the nickel perovskites, and 
in a broader perspective, in correlated electron systems. 



II. EXPERIMENTAL 

The pressure dependent XAS measurements at the Ni 
K-edge ( 8345 eV) were performed at the dispersive XAS 
beamlin^iSlin! of the LNLS (Laboratorio Nacional de Luz 
Si'ncrotron, Campinas, Brazil). A fine grained high qual- 
ity YNiOs polycrystalline powder sampl^ was loaded in 
the 125 fim diameter hole of an iconel gasket mounted on 
2.4 mm thick diamond anvils, with a cullet of 300 fim. 
The spectra were measured using a Si(lll) bent crys- 
tal monochromator that was selecting a band-pass of 
about 800 eV from the white beam and focusing it into 
a Gaussian spot of 150 /xm -FWHM at the sample po- 
sition. The gasket cut the tails of the Gaussian beam. 
However, the beam position for each energy was laterally 
dispersed in less than 30 /im, such that the full energy 
band-pass were transmitted through the gasket with al- 
most the same photon flux. The reference flux, Iq, was 
measured through a flat piece of glass in order to sim- 
ulate the average attenuation without introducing new 
features in the spectra. At each pressure, the cell was 
realigned at the optical focus. In dispersive XAS there is 
no mechanical movement of the optics during data collec- 
tion and the whole spectrum is observed at once. It then 
is possible to screen the cell orientations in a relatively 
short time, to find an optimal position. The remaining 
glitches were then removed by small rotations of the cell 
with respect to the polychromatic X-ray beam. Using 



a gas membrane-driven mechanism, the pressure was in- 
creased by steps of about 2 GPa, measured using the 
ruby methocpi'. Pressures up to 20 GPa were applied 
using silicone oil as pressure medium. Above 15 GPa the 
non-hydrostatic components of this medium may lead to 
pressure deviations up to 10% over a diameter of 150 fim. 
The precision on the pressure is then around 0.5 GPa up 
to 10 GPa. Above this value, an error bar of 10%, associ- 
ated to the non-hydrostaticity, is estimated for the abso- 
lute pressure scale. Since XAS probes the effective aver- 
age contribution from all grains randomly oriented over 
the pressure gradient in the sample area and due to the 
order of magnitude of the effects (of the order of 10~^), 
XAS is not very sensitive to a limited non-hydrostaticity. 

EXAFS (Extended X-ray Absorption Fine Structure) 
data free of Bragg peaks were collected up to about 
10. 5A . This data range limits the minimum differ- 
ence between two close distances that can be resolved 
in the EXAFS analysis to the AR ~ 0.15 A. In the 
YNiOs under ambient conditions the average Ni-0 dis- 
tance is 1.994 A and the separation between the largest 
and shortest distance is only 0.039 A^. The restricted 
EXAFS k-range (~10 A~^) characterizes then a low- 
resolution study in R-space. The Ni-0 bond lengths are 
not distinguishable and bond length differences within 
the two Ni sites appears as a static disorder contribu- 
tion to the total bond length dispersioiP^. The EXAFS 
analysis gives only the average distance. The limited k- 
range broadens the Ni-0 peak but does not change its 
average position. Shifts in this average distance with in- 
creasing pressure can be tracked precisely, provided that 
the analysis procedure is strictly the same for the whole 
series of measurements. The data were analyzed using 
the Athena/ Artemis packagd^l, xhe EXAFS oscillations 
were extracted following the standard procedure of the 
Athena code. The signal corresponding to the oxygen 
coordination shell was selected by Fourier filtering and 
the structural parameter were deduced from fitting pro- 
cedures. In the fitting procedure the number of free pa- 
rameters was limited by the useful range and the interval 
corresponding to the selected signal in the real space. Af- 
ter a first screening of the individual data sets, all data 
were fitted together using the option of multiple set anal- 
ysis. This option increases the number of independent 
points over the number of free parameters. By setting 
the same (fitted) origin to the k-scale, it reduces sub- 
stantially the uncertainty on the relative decrease of the 
Ni-0 distances with pressure. The number of neighbors 
and amplitude factor were fixed. The only amplitude pa- 
rameter left for each data set is the Debye- Waller term tr, 
that describes the dispersion of the Ni-0 bond lengths. 
The changes in Ni-0 bond length SR^i-o and the bond 
length dispersioncT are obtained with error bars around 
0.005A and 0.02A, respectively. 

For the XANES (X-ray Absorption Near Edge Structure) 
range the normalization procedure is the same for all 
pressures. In order to avoid spectra deformation we limit 



the data handling to simple operations: normalization 
consists essentially in the subtraction of a straight line 
for background and to a scaling of the edge jump to 1 in 
the range 150-250 eV above the edge. Variations in the 
straight line subtraction and the exact position of the 
normalization zone affect to a small extent the overall 
shape of the spectra but very little the relative intensity 
of the structures. Due to the monochromator bandwidth, 
detector spatial resolution and the core-hole lifetime, the 
experimental resolution for the XANES experiments is 
around 1.5 eV. However, a measure of the shift of the 
edge energy with pressure is not limited by the experi- 
mental resolution but by the stability of the spectrom- 
eter and counting statistics. In a dispersive XAS setup 
the monochromator does not move allowing a very high 
energy stability. During the whole experiment, a sharp 
Bragg peak outside the analysis range allowed us to ver- 
ify that the energy stability was better than 50 meV, in 
accordance with previous measurement^2ni_ 

The XANES features were compared to ab initio full mul- 
tiple scattering calculations using the FDMNES codtP 
for Ni-centered clusters and with ato mic p ositions given 
by reported crystallographic structures!'^. 



III. RESULTS 



Figure 1-a shows the absorption spectrum at ambient 
pressure and the EXAFS signal for selected pressures (in- 
set). Figure 1-b gives the modulus of the Fourier trans- 
form (FTM) of these signals. This representation gives 
the pseudo-radial distribution function (not corrected for 
EXAFS phase shifts) around the average Ni atom. The 
position of the most prominent peak around 1.5 A corre- 
sponds to the average Ni-0 bond length. The next two 
peaks correspond essentially to the Ni-Y bondings and 
to the Ni-Ni backscattering and Ni-O-Ni multiple scat- 
tering, respectively. These contributions reveal the sensi- 
tivity of EXAFS to the next-nearest neighbors geometry 
and, so, to the bond angles among octahedra. 

The progressive shift of the first peak of the FTM to- 
wards shorter distances (Fig. 1-b) yields the contraction 
AR in the average Ni-0 bond lengths. Figure 2 shows the 
pressure dependence of the Ni-0 distance contraction, 
obtained from the quantitative analysis after selection of 
that peak. Between 2 and 13 GPa, AR decreases al- 
most linearly with the applied pressure. The contraction 
AR up to 13 GPa gives w 0.033 A. Such a contraction 
represents a NiOq relative volume decrease of 5.0%. In 
the same range the unit cell volume measured by diffrac- 
tion also drops by 5.0%^^. The monotonic decrease of 
the Ni-0 average distance stops at around 13 GPa. In 
a transient pressure range (13 GPa < P < 16 GPa), the 
average Ni-0 bond length and the distance dispersion 
remain almost unchanged. Above 16 GPa, average Ni- 
O bond length resumes shrinking. The total contraction 
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Figure 1: a: typical Ni K-edge XAS spectrum of YNiOz- In- 
set: EXAFS signal in the available k-range for different pres- 
sures, b Fourier transform modulus of the EXAFS oscillations 
for selected pressures. 



up to 20 GPa is w 0.043 A, or 6.6% in terms of volume 
decrease. The increase of the amplitude of the first peak 
(Fig. 1-b) corresponds to a decrease in the Ni-0 bond 
length dispersion. The largest increase occurs from 7 
to 13 GPa. Within this range, it the average iVzOg units 
tend toward a partial symmetrization. The evolution of 
the Debye-Waller term (Fig. 2, inset) shows that indeed 
the largest decrease in the bond dispersion occurs from 7 
to 13 GPa, confirming this qualitative outcome. 

The strongest modifications in the medium-range struc- 
tures take place within the pressure range 13-16 GPa. 
As seen in figure 1-b, the intensity of the Ni-Y peak in- 
creases while those of the Ni-Ni and Ni-O-Ni peaks de- 
crease. These modifications point out inter-octahedra re- 
arrangements. Qualitatively, the evolution of position of 
the Ni-Y peak towards low R is easily interpreted in terms 
of shortening of the Ni-Y distances. The "shift back" to- 
wards larger distances observed at pressures higher than 
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Figure 2: Pressure dependence of the average Ni-O bond 
length contraction Ai? obtained from the EXAFS analysis. 
The line is just a guide for the eyes. Vertical lines separate 
different regions discussed in the text. Inset : o"'^ as a function 
of the pressure. 



Figure 3: XANES spectra in YNiOa for selected representa- 
tive pressures. A is the pre-edge, B the main edge structure 
and C the shoulder associated with middle range effects (see 
text). Lower inset : zoom on the C feature. Upper inset: 
maximum of the derivative at the edge as a function of the 
pressure. 



15 GPa does not fit with such a simple scheme. One 
may note that aheady at 15 GPa the peak is substan- 
tially enlarged on the high R side. As this enlargement 
occurs at the pressure for which the monoclinic to or- 
thorhombic transition is expected, it could be a finger- 
print of the phase transition, and consequently of some 
inter-octahedra rearrangement. 

Figure 3 shows the NiK-edge XANES spectra of YNi03 
at selected pressures. The pre-edge feature A, which 
probes the Ni 3d states hybridized with Ni 4p states, 
does not change within the experimental precision. Two 
main changes are observed when the applied pressure is 
increased: a shift of the absorption threshold and subtle 
modifications of the spectral features B and C. 

Edge shifts are primarily associated to changes in the 
formal valence. However, as here the Ni keeps the form 
Ni'^'^ , these shifts have essentially a structural origin. As 
coordination bond lengths decrease, the band energy will 
increase as a result of the higher overlap of the electrons 
density of neighboring atoms. This increased overlap is 
less significant for localized Is levels than for 4p delocal- 
ized ones and, consequently, the absorption edge shifts 
towards higher energies. For small shifts the bond length 
reduction (SR) and t he ass ociated edge shift (SE) are 
almost linearly relatecP^^. However, due to the mix- 
ture of close Ni bonds such linear relationship does not 
strictly apply in YNiO^, and the evolution of the edge 
yields only a qualitative analysis. 

The edge shift SE is measured by the position of the 
derivative maximum. The observed shift rate, 6E/6P, 
as a function of applied pressure is not constant over 
the whole pressure range (Fig. 3- inset). Up to 13 GPa, 
SE varies linearly with the applied pressure and SE/SP 
is about 60 meV/GPa. Within the range 13 GPa < P 
< 16 GPa, the edge position position shifts at a much 



slower rate. A new increase in SE/SP above 16 GPa 
indicates the recovering of a regime of bond contraction. 
All these outcomes are in full qualitative agreement with 
the EXAFS analysis. 

Within the range 13 GPa < P < 16 GPa, a subtle raise 
of a shoulder C about 10 eV above the main structure 
B is observed. The emergence of simi lar shoulder was 
reported in NdNiOz and PrNiO^^^. In these studies 
it has been clearly associated to thermal induced met- 
allization. These authors relate the appearance of the 
shoulder C to changes in the local atomic environment 
around Ni at the IM transition. The physical origin of 
this spectral feature will be discussed below. Neverthe- 
less, we can infer that in the present case the shoulder G 
must also be associated to metallization. We note that its 
appearance takes place in a narrow range around 14 GPa, 
where metallization is indeed expectecP^, but also where 
the EXAFS analysis shows that there is no significant 
change in the Ni-0 bond lengths. 

We simulated the Ni K edge XANES spectra for YNiOz 
using ah initio full- multiple scattering calculationJ^with 
atomic positions given by crystallographic structures re- 
ported by neutron diffractiorpJ^. We observed that the 
simulations using the monoclinic or the orthorhombic 
crystallographic structures for YNiOs lead to identical 
spectra, and do not show the expected C shoulder (Fig. 4). 
Similar comparison of the calculated spectra in PrNiOs, 
for the reported orthorhombic structure in the metallic 
and insulator phases, led Acosta and coworkerJ^^ to de- 
duce that the changes in the local atomic structure are 
not reflected in the average crystalline structure. Our re- 
sults clearly support this view. In addition, the absence 
of the structure G in our calculation based on the mon- 
oclinic structure shows that this feature is not related 
to the existence of two Ni sites. This suggests that G 
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Figure 4: ab initio simulations of the XANES spectra for a Ni 
centered 61 atoms-cluster of YNiO^ in the monoclinic (plain 
line), orthorhombic (dashed line) and LaNiO^-Wke rhombohe- 
dral (dotted line) structures Inset : simulations for the rhom- 
bohedral structure and different cluster sizes . 

comes from some characteristics of the local organization 
beyond the coordination shell. 

To gather elements about the origin of this feature, we 
went further simulating the XANES spectra based on 
a rhombohedral structure. Even if not reported for 
YNiOs, rhombohedral structure corresponds to the most 
symmetric structure in the RNiO^ series. Importantly, 
in this structure the tilt angles among adjacent octahedra 
are quite different compared to those found in monoclinic 
and orthorhombic YNiO^. This provides a track to in- 
terpret our data. The cluster was built from cell param- 
eters deduced from those of LaNiO^ substituting La by 
Y in the cell. In order to facilitate the comparison with 
the experimental features in YNiOs, we use an isometric 
expansion of these parameters keeping the cell volume 
that of the YNiO^. In this rombohedral structure the 
nickel atoms have a unique regular octahedral site and 
the superexchange Ni-O-Ni angle is 165 degrees. As for 
all previous calculations, the cluster contained 61 atoms. 
The shoulder C is present around 10 eV above B (Fig. 4). 
We then performed simulations using the same LaNiO^- 
rhombohedral structure but with decreasing cluster size 
(Fig. 4 inset). We checked that the C shoulder does not 
correspond to a specific path of multiple scattering within 
the octahedron. This structure emerges when the cluster 
considered in the simulations includes at least 46 atoms 
(cluster size ss 4.8 A), i.e. half of the neighboring oxy- 
gen's octahedra. 



IV. DISCUSSION AND CONCLUSION 

We used X-ray absorption spectroscopy to describe the 
evolution of the local and middle range order under ap- 
plied pressure. A first regime of monotonic contrac- 



tion and symmetrization of the average NiOe units takes 
place up to 13 GPa. It stops in a limited pressure do- 
main around 14 GPa, before resuming above 16 GPa. 
In this narrow pressure range, crystallographic modifica- 
tions basically occur in the medium/long range, not in 
the NiOe octahedron, whereas the evolution of the near- 
edge XAS features can be associated to metallization. 
One important outcome in our experiment is the emer- 
gence of the shoulder C in the XANES around 14 GPa. 
Such emergence is reproduced in the simulations when 
the cluster radius reaches a critical value of « 4.8 A. It 
is then not related to the actual very local geometry or 
symmetrization of the Ni sites but to middle range ef- 
fect. This view is supported by the experimental out- 
come that the appearance of the shoulder C occurs in a 
narrow pressure range around 14 GPa where no signifi- 
cant modifications in the average NiOe octahedron are 
observed, whereas significant modifications seem to take 
place in the middle range order. We should point out 
that in t he exp erimental studies reported in PrNiOs and 
NdNiO^^, the shoulder C rising at the IM transition 
is more intense than in our study. We attributed that 
to our limited resolution, where a compromise had to be 
found between good resolution for XANES and energy 
range for EXAFS. As a consequence, the absolute inten- 
sity of the experimental features C is not well reproduced 
in the simulations. 

The main effect of pressure is to decrease bond lengths 
and straighten bond angles leading to an increased band- 
width (W) and to a correlated decrease in TnJ^^. Tim 
goes down to room temperature at 14 GPa owing to 
the increased bandwidth, which changes the ratio W/U, 
where U is the on-site Coulomb energy. U is not expected 
to significantly modify around 14 GPa because it con- 
cerns much more localized 3d orbitals and, in addition, 
the bond length modification is very small. Although 
electron-electron correlation (U) represents the main en- 
ergy in these strong correlated systems, the increase in 
the band width (W) seems to control the electronic tran- 
sition. 

The presence or not of the shoulder C around 15 eV 
above the absorption threshold in the XANES spectra 
of rare earth nickelate turns out to be related to the abil- 
ity of the Ni 3d eg and O 2p orbitals to overlap to form a 
band. Besides the Ni-0 bonding interaction this ability is 
strongly related to the Ni-O-Ni angle. The Ni-O-Ni angle 
crossover for the IM transition at room temp erature and 
pressure in RNiO^ is around 156 degreea^ESl. In both the 
monoclinic and orthorhombic crystallographic structures 
given by neutron diffraction the Ni-O-Ni angle is around 
146 degrees in YNiO^, while this angle is aroimd 157 
degrees in PrNiO^ and NdNiO^. The XANES results 
are in perfect agreement with the EXAFS outcome that 
around 14 GPa the reduction of the distance within the 
octahedra slows down, while a significant reorganization 
takes place for the Y and Ni next neighbors. The drop in 
the lattice parameter c identified by Garcia-Munoz and 
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coworkers at 14 GPcP^ corresponds more likely to inter- 
octahedra rearrangements than to a modification in the 
local geometry of the NiO^ units. 

In conclusion, pressure dependent X-ray absorption spec- 
troscopy in YNiOz shows that, after an initial shrinking, 
there are no significant modifications in the average NiOq 
octahedra around 14 GPa, where a signature of electronic 
delocalization is observed. This experimental outcome, 
supported by ah initio XAS calculations, provides evi- 
dences of the bandwidth driven nature of the pressure 
induced IM transition in the same way that it has been 
recently reported in LaMnO^^. The electronic delocal- 



ization is essentially synchronized with the opening of the 
Ni-O-Ni angle, an not to a sudden modifications in the 
local geometry of the NiO^ octahedra. 
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